Effective Hamiltonians in quantum optics: a systematic approach 
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We discuss a general and systematic method for obtaining effective Hamiltonians that describe 
different nonlinear optical processes. The method exploits the existence of a nonlinear deformation 
of the usual su(2) algebra that arises as the dynamical symmetry of the original model. When some 
physical parameter, dictated by the process under consideration, becomes small, we immediately get 
a diagonal effective Hamiltonian that correctly represents the dynamics for arbitrary states and long 
times. We extend the technique to su(3) and su(N), finding the corresponding effective Hamiltonians 
when some resonance conditions are fulfilled. 
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^0 ' Quantum optics provides the ideal arena to deal with the interaction of radiation and matter. Indeed, by using 
\ standard techniques one easily gets a system of nonlinear coupled equations that govern, under the approximation 
of negligible damping, the interaction of a discrete set of field modes with an ensemble of atoms ^ ^ Q. This 
is a very general description, but too involved to draw any immediate physical conclusion. Nevertheless, as soon as 
»^ , one realizes the existence of fast and slow variables with a large difference in time scales, it is rather obvious that a 
simplified formulation can be derived from first principles, at least for some limiting cases. 

In most nonlinear optical phenomena the fields and the atomic transitions are usually far from resonance 
• In fact, no appreciable population redistribution is brought about the irradiation even of intense fields, and only 
\ minute fractions of light intensities are absorbed. In the reference frame of the external field oscillations, the atomic 
^NJ ' polarization is then a fast variable controlled by the slow motion of the electromagnetic field amplitudes and follows 
' its evolution adiabatically. Under this assumption we can eliminate the atomic degrees of freedom and are left with a 
: small number of equations for the fields alone. Due to the elimination procedure, these equations appear as suffering 
^ ^ a nonlinear field-field interaction and can consequently be reinterpreted as the Heisenberg equations of motion for the 
^ ' , , field operators under the dynamics of an effective nonlinear Hamiltonian. In other words, such an effective Hamiltonian 
' acts only on the slow degrees of freedom and, because it correctly describes the slow motion, it incorporates the effect 
d of the coupling of these slow degrees of freedom with the fast ones. 

^ However, this adiabatic elimination (that was first used to study the parametric oscillator presents some draw- 

backs. First, it does not provide a general prescription for finding effective Hamiltonians, since the particular details 
strongly depend on the model considered. Second, it could become very cumbersome, and explicit but enormously 
complicated expressions for the different orders of approximation can be found in many original publications ^. 
^ . , Third, and even worse, the procedure is not uniquely defined: depending on the term eliminated, the outcome of the 
■ final Hamiltonian could be different |l2|. 

For these reasons, other methods of deriving effective Hamiltonians exist. For the field of quantum optics, the ones 
devised in Refs. jl^ and are especially germane. A thorough review of the different ways of formally constructing 
effective Hamiltonians may be found in Ref. p^ . Roughly speaking, all the methods have in common that at some 
point in their implementation one applies a unitary (or canonical) transformation to the total Hamiltonian and keeps 
only terms up to some fixed order. It is worth emphasizing that this technique is standard in condensed matter 
physics fl^ and important examples of its application may be given fl^ . 

In consequence, it seems pertinent to find a setting to support this usual approach of effective Hamiltonians: it is 
the main goal of this paper to provide such a setting by resorting to some elementary notions of group theory. To 
this end, we first note that most of the effective Hamiltonians in quantum optics contain cubic or higher terms in 
creation and annihilation operators. Among others, typical examples are fcth harmonic generation, fc-wave mixing, and 
generalized Dicke models. The key point for our purposes is the recent observation that the common mathematical 
structure underlying all these cases is a nonlinear or polynomial deformation of su(2), which arises as the dynamical 
symmetry algebra of the corresponding Hamiltonian. This nonlinear algebra has recently found an important place 
in quantum optics ^ ^ because it allows us to handle the problems in very close analogy with the usual 

treatment of an angular momentum. In particular, we get a decomposition of the Hilbert space into direct sums of 
invariant subspaces and the dynamical problem generated by the corresponding Hamiltonian can be reduced to the 
diagonalization of a finite-dimensional matrix. 
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This is a considerable achievement by itself, but unfortunately it is im pos sible to obtain analytic expressions for 
the eigenvalues and eigenstates of those matrices. For this reason, in Ref. we have devised a Lie-like method |2j] 
that allows one to get approximately effective Hamiltonians that can be diagonalized in an exact form. The idea is 
to apply a small nonlinear "rotation" to the original Hamiltonian: although, in general, its action is rather involved, 
when some physical parameter (dictated by the particular model under consideration) becomes small, it generates an 
effective Hamiltonian that is diagonal and represents correctly the dynamics for arbitrary states and even for long 
times. The method appears then as the natural and systematic tool for finding effective Hamiltonians, provided the 
existence of this nonlinear su(2) dynamical algebra. 

Furthermore, it is worth noting that the su(3) algebra is the natural extension of su(2) to study the dynamical 
evolution of three-level systems. Far from being an exotic curiosity, this evolution is central to the discussion of many 



fascinating problems, such as two-photon coherence |23, 26 , resonant Raman scattering G8l, superradiance |2 



and three-level echoes [^0|, |3l|]. As one would expect, when these three-level systems interact with quantum fields, a 
nonlinear deformation of su(3) naturally emerges. Much in the same way, a su(N) deformed structure arises when 
considering the interaction of A''-level systems with quantized fields. 

In this paper we apply our method to some nontrivial examples: a collection of two-level atoms dispersively 
interacting with a quantum field in a cavity; and S and A configurations of three-level systems interacting with fields 
under different resonance conditions. We also obtain effective Hamiltonians for multilevel systems under A^-photon 
resonance conditions. Finally, wc present how our method could work for more general situations, which would greatly 
facilitate the physical applications of this technique to many quantum optical problems. 

II. EFFECTIVE HAMILTONIANS FOR NONLINEAR SU(2) DYNAMICS 

A. Motivation of the method 

To keep the discussion as self-contained as possible and to introduce the physical ideas underlying the method, let 
us start with the very simple example of a particle of spin j in a magnetic field. The Hamiltonian for this system has 
the following form (in units h = 1, which will be used throughout all this paper) 

H^iuS3+g{S+ + S-), (2.1) 

where g is the coupling constant and the operators 5*3,5'+, and 5_ constitute a {2j + l)-dimcnsional representation 
of the su(2) algebra, obeying the usual commutation relations 

[53,5±]=±5±, [5+,5_]= 253. (2.2) 

In the traditional angular momentum basis \j, m) (m = — j, — j -I- 1, ... ,j — 1, 7) the operator 53 is diagonal 

'5'3|j, = "i), (2.3) 

while the action of the ladder operators S± is given by 



S±\i,m) = ^/{iTm){j ±m+l)\j,m±l). (2.4) 

The Hamiltonian (|]^) belongs to the class of the so-called linear Hamiltonians |Q and admits an exact solution. A 
convenient way of working out the solution is to apply the rotation 

[/ = exp[a(5+ - 5_)] , (2.5) 

and recalling that e^Be^"^ = B + [A^ B] + ^[^, [A, B\\ + . . ., the rotated Hamiltonian, which is unitarily equivalent 
to the original one, becomes 

H = UHU'^ = [wcos(2a) + 2gsin(2a)]53 + i[2.gcos(2a) - w sin(2a)](5+ + 5_). (2.6) 

Now, the central idea is to choose the parameter a so as to cancel the nondiagonal terms appearing in Eq. (^.(j). 
This can be accomplished by taking 

tan(2a) = — , (2.7) 

UJ 
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and the transformed Hamiltonian reduces then to 



^^off = w\/ 1 + ^ S3. 



(2., 



Since this effective Hamiltonian is diagonal in the angular momentum basis, the dynamical problem is compl etely 
solved. The crucial observation is that when g <^ m we can approximate Eq. (2.7) by a ~ g/i^, and then Eq. ( |2.5| ) 
can be substituted by 



U ~ exp 



-iS+-S-) 



This small rotation approximately (i.e., up to second-order terms in g/cu) diagonalizes the Hamiltonian ( p.l[ ), giving 
rise to the effective Hamiltonian 



(2.9) 



2^ ) ^3 

UJ 



(2.10) 



which obviously coincides with the exact solution (|2.8D after ex pand ing in a series of g^ juP'. A direct application of 
the standard time-independent p ertur bation theorypsf to Eq. ( ^.l[ ) leads immediately to the same eigenvalues and 
eigenstates that the Hamiltonian ( 2.10 ) gives in the same order of approximation. However, we stress that our method 
is fully operatorial and avoids the tedious work of computing the successive corrections as sums over all the accessible 
states. 



B. Nonlinear small rotations 



Having in mind the previous simple example, let us go one step further by treating the more general case of a 
system that admits some integrals of motion Nj and whose interaction Hamiltonian can be written as 

i?i„t - A + g{X+ + X_), (2.11) 

where (7 is a coupling constant and A is a parameter usually representing the detuning between frequencies of different 
subsystems (although it is not necessary). The operators X± and X3 maintain the first commutation relation of su(2), 
[X3,X-|-] = ±X±, but the second one is modified in the following way 

[X+,X_] = P(X3), (2.12) 

where PiX-^) is an arbitrary polynomial function of the diagonal operator X^ with coefficients perhaps depending 
on the integrals of motion TV^ . These commutation relations correspond to the so-called polynomial deformation of 
su(2). Such nonlinear algebras were discovered by Sklyanin jsj] and Higgs ]35[| and have already played an important 
role in quantum mechanics |3^, |3^, |3^ . 

Now suppose that for some physical reasons (depending on the particular model under consideration) the condition 

£ = I « 1 (2.13) 



is fulfilled. Then, it is clear that (2.11) is almost diagonal in the basis that diagonalizes X3. In fact, a standard 
perturbation analysis immediately shows that the first-order corrections introduced by the nondiagonal part g{X^ + 
X_) to the eigenvalues of X3 v anish and those of second order are proportional to e. Thus, we apply the following 



unitary transformation to (2.11) (which, in fact, is a small nonlinear rotation) 

[/ = cxp[e(X+ -X_)]. (2.14) 

After some calculations we get, up to order e'^, the effective Hamiltonian we are looking for 

2 

i/eff = A X3 -f ^PiXs), (2.15) 



which is diagonal in the basis of eigenstates of X3. Then, the evolution (as well as the spectral) problem is completely 
solved in this approximation. Besides the advantage of having the effective Hamiltonian expressed in an operatorial 



4 



form, the method has the virtue of generahty, since it is vaUd for any model whose Hamiltonian could be written 
down in terms of the generators of a polynomial deformation of su(2). 



Our technique also provides a valuable tool for obtaining corrections to the eigenstates of (2.11). Indeed, it is easy 
to realize that these eigenstates can be approximated by 



(2.16) 



where \m) denotes an eigenstate of and U is the corresponding small rotation. Since U and |m) do not depend 
on time, the operator U can be applied to |m) as an expansion in e. For example, the eigenstate |5'„i) up to order 
takes on the form 



1-e (X+-X_) 



(1 + 2X+X^ 



X' 



(2.17) 



This representation is especially advantageous when the state space of the model is a representation space of the 
deformed su(2) algebra that is constructed in the usual way by the action of the raising operator X.). ; i.e., jm) ~ X™|0), 
where |0) is a lowest weight vector fulfilling the standard condition X_|0) = 0. 

This general procedure shows that can always adiabatically eliminate all nonresonant transitions and work with 
an effective Hamiltonian containing only (quasi) resonant transitions. The effect of nonresonant terms reduces to a 
dynamical Stark shift (which can have a quite complicated form). Obviously, transformations generating effective 
Hamiltonians also change eigenfunctions, but the corrections are of order e and do not depend on time. Such 
corrections correspond to low-amplitude transitions that take place in the case of nonresonant interactions. We shall 
elaborate on these topics in the next Sections. 



C. An example: dispersive limit of the Dicke model 



As a relevant example, let us apply our method to the well-known Dicke model that describes the interaction of a 
single- mode field of frequency Wf with a collection of A identical two- level atoms with transition frequency wq- Making 
the standard dipole and rotating- wave approximations, the model Hamiltonian reads as [Eol 



H = ujio) a + LOaS^ + g{aS+ + a^5_). 



(2.18) 



where a(a^) are the annihilation (creation) operators for the field mode, and {S±, S3) are collective atomic operators 
forming an {A + l)-dimensional representation of su(2). 

By introducing the excitation number N — a^a + S3, which is an integral of motion, we can recast ( 2.18 ) as 
H ^Ho+ Hint with 



Ho = WfiV, 
Hint = A 53 + giaS+ + a'^S^), 



(2.19) 



the detuning being A — luq — ujf. 

We assume now that the dispersive limit holds l41[; i.e., 



|A| > AgVf 



1, 



(2.20) 



where n is the average number of photons in the field. If, after our previous discussion, we introduce the deformed 
su(2) operators as 

X+^aS+, X_=at5_, X3 ^ S3, (2.21) 
then we immediately get that the small nonlinear rotation (2.14) transforms the interaction Hamiltonian ( 2.19| ) into 



H,s = AS3 + ^[Si - 2(ata + 1)^3 _ C2], 



(2.22) 



where C2 — A/2{A/2 + 1) is the value of the Casimir operator for su(2). The effective Hamiltonian (2.22) was 
previously obtained in Ref. by quite a different method (see also Ref. ^^) and, due to the presence of the 
nonlinear term 6*3 , has been considered as a candidate for the generation of squeezed atomic states H] . 
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III. EFFECTIVE HAMILTONIANS FOR NONLINEAR SU(3) DYNAMICS 
A. Three-level systems interacting with quantum fields 

The method of small rotations discussed in the previous Section can be applied not only to Hamiltonians having 
a nonlinear su(2) structure, but also to more complicated systems. In this Section we shall focus on Hamiltonians 
that can be represented in terms of the su(3) algebra. This structure naturally arises when dealing with systems with 
three relevant levels. It is well known that in this case three possible configurations (commonly called S, V, and A) 



are admissible (see Ref. |45 for details). 

For definiteness, we shall consider the interaction of A identical three-level systems in a cascade or S configuration 
(i.e., with associated energies Ei < E2 < E3 and allowed dipole transitions 1^2 and 2^3, but not the 1^3) 
interacting with a single-mode quantum field of frequency uj{. The Hamiltonian of this model is 

He — ^^fiold + -f^atom + ^^int, (3.1) 

with 

Wfieid = ujfa^a, 

Hatora — EiS^^ + i?2<5'^^ -|- E^S^^ , (3.2) 

As usual, the three diagonal observables 5*' — {i — 1,2,3) measure the population of the ith energy level, 

while the off-diagonal polarizations S'*-' = \j){i\ generate transitions from level i to j. They satisfy the commutation 
relations 

[S'',S'''] = S,kS''-5aS''', (3.3) 

which turn out to be those of the u(3) algebra. In fact, they form a {A + 1){A + 2)/2-dimensional representation of 
the u(3) algebra. 

Because the sum S^^ + + S^'^ = A is an obvious integral of motion that determines the total number of atoms, 
only two of the populations can vary independently. For this reason it is customary to introduce the two traceless 
operators 

^12 = 1(^22 _ ^11)^ = 1(^33-522), (3.4) 

that represent population inversion between the corresponding levels. Furthermore, to emphasize the idea of transition 
between levels, it is usual to define the operators S^f — S*'-* and S'j' — 5^* (for i < j). Then the eight independent 
operators {S^^,Sl^) satisfy the commutation relations of su(3). Moreover, {S}^,Sl'^) and {S'^,S^'^) form two su(2) 
subalgebras (each one describing the transitions 1 <-» 2 and 2^3). Nevertheless, transitions 1^2 and 2 3 are 
not physically independent, since we have 

[Sl\Sf]^-Sl\ [Sl\Sl']^S]_\ [Sl^,S^_^]^0. (3.5) 

It seems natural also to introduce the deformed su(3) algebra as 

j^ll _ ^11 j^22 _ g22 j^33 _ g33 



(3.6) 



Then the Hamiltonian (3T) can be recast as Hs — Hq + ifint, with 



Ho = uJiNs, 



(3.7) 



i?i„t = -Ai2 X" + A23 + gi2(X|2 + Xl2) + g23(Xf + Xl^), 
where we have used the conserved excitation number iVs = a^a + S^^ — S^^ and we have introduced the detunings as 



A12 — E2 — El — uji, A23 — E3 — E2 — 



(3.8) 
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The operators X"^^ satisfy the usual su(3) commutation relations, with ( |3.5D recast as 

[X'^,Xf] = -Yl\ [Xl2,X23] ^y^s, [X|2,X23]=0. (3.9) 

where 

Yl^ - a'S'^. (3.10) 

However, we have to modify some of them in the following way: 

[X'l ,X'1]^P{X'\ X" ) , ■. ^- \^Q{X'\X" ), (3.11) 

where P(X",X'^'"') and Q{X" , X'^'') are polynomials of the diagonal operators X" {i = 1,2,3) and define, then, a 
polynomial deformation of su(3). 



1. Effect of a far-off resonant level 

The dynamics generated by su(3) is obviously richer than that of su(2), since a greater number of physical degrees 
of freedom are now available. To see how our method works in this case, let us assume that one of the transitions, 
say the 2 <-!■ 3, is (quasi) resonant with the field; i.e., 

IA23I <^523V^^, (3.12) 

while the transition 1 ^ 2 is far-off resonant 

IA12I > A9i2Vn+l. (3.13) 
It is clear from our previous analysis that the small nonlinear rotation 

[/12 -exp[ei2(X|2-Xl2)], (3.14) 

with £12 = gi2/Ai2 <C 1, eliminates the interaction term gi2{X^ + Xl^), representing the nonresonant transition 
1 <-» 2. In such a way, we obtain the effective Hamiltonian 

i/eff - A12 X" + A23 + 523(Xf + - 512523 (^^3 ^ ^^3) 

^12 

+ |lLp(x",x22) + ^H5Hi([xi^x23] + [Xf,Xl2]). (3.15) 

Ai2 Ai2 

The remarkable point is that by eliminating the transition 1 <-> 2, we have generated an effective transition 1 <-> 3 
(represented by the operators Y^^), which was absent in the initial Hamiltonian. Nevertheless, this transition is also 
nonresonant due to conditions ( |3.12 ) and (3.13) and, accordingly, can be eliminated by the following transformation 



C/i3 = exp[£i3(r^3 - ri^)], (3.16) 
where the parameter £13 in the above equation must fulfill 

512523 ^ , CQ 17\ 

= A,2(A,2+A23) « ' ^'-''^ 

In this particular case the polynomial function P{X'", X^^) is 

P(X",X22) = Sl^Sl^ + a^a{S^^ - 5"). (3.18) 

If initially level 1 is unpopulated we have = (which will be conserved, since there are no transitions to the level 
1), and 

F(X",X22) = 522(ata + i). (3.ig) 
In consequence, the effective Hamiltonian, taking into account the existence of a far-off resonant level, has the form 

ffeff = A23 S'^+g23iaSf+a^Sl^) + + (3.20) 

^12 

This means that the far-lying level produces a mark in the system in the form of a dynamical Stark-shift term. It 
is interesting to observe that due to this Stark shift, the initially nonresonant transition 2 <-> 3 becomes resonant in 
a subspace with some fixed photon number. This opens the possibility of separating n-photon field states from an 
initial coherent state interacting with a collection of three-level atoms just by choosing a suitable relation between 
the detunings and the interaction constants, and projecting to the second level at appropriate moments. 
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2. Two-photon resonance 

Let us now envisage the different situation in which the two-photon resonance condition between levels 1 and 3 is 
fulfilled; i.e., E^ — Ei = 2wf . This means that A12 = — A23 and the transition generated by the operators Y^^ cannot 
be removed. However, the term g23{X^ + X'^), which gener ates (nonresonant) transitions between levels 2 and 3, 



can be eliminated by a transformation U23 analogous to ( |3.14 ) with £23 <C 1. The transformed Hamiltonian becomes 
then 

iJeff = -A12 {X'^+X'')-^iYl' + Yl') 

A12 

+ |kp(X",X22)_|iip(X22,X33), (3.21) 
A12 A12 



where P{X-^^ , X^'^) and P{X^^ , X^^) are defined according to Eq. ( 3.18| ). Finally, if we further impose the absence of 
initial population in level 2, we obtain 

H,^ = »(a25i3 + at'5i3) 
A12 

+ iSl' + A/2) [igl,/Au - 5?2/Ai2)a^a + 523/A12] + A^a^a, (3.22) 

^12 

which is the effective two-photon Dicke Hamiltonian including the dynamical Stark shift obtained by Puri and Bul- 
lough (see also Rcf.|l|). 

B. The dispersive limit of the A configuration 

Let us consider for a moment the case of the A configuration, in which the allowed dipole transitions are now 1 <-> 3 



and 2 <-!■ 3, but not the 1 <-> 2. The Hamiltonian governing the evolution is still of the form (3.1), but now with 

Hint = gisiaSl^ + a^Sl^) + 523 (a^f + a^S^_^). (3.23) 
By using the integral of motion N_\ ^ a^a + S^^ we can rewrite 

Hint = ~A3iX" - A32X22 + (Xf + X'J) + 323 (^f + Xl^), (3.24) 



where, as in Eq. (3^), we have introduced the deformed su(3) operators as 



Xl^ = aSl^, Xf = aSf. (3.25) 



These deformed generators satisfy a set of commutation relations similar to (3.9), but instead of ( 3.10| ) we must use 



yi2 ^ (^33 _ ^tfl), (3.26) 

which is the mathematical reason for the well-known different behaviours exhibited by S and A configurations. 
Let us focus on the dispersive regime, when 



IA13I > AgisVn + T, IA23I > A923\/n+ 1. (3.27) 
Then a couple of small rotations eliminate the far-off resonant transitions 1 3 and 2 <-> 3, obtaining 
Hcs — — A3iS'^^ — A325^^ 

+ |ii[(5ll + 1)^33 ^ ^ta(533 _ ^11)] + |k[(^22 ^ ^)^33 _^ ^1^(^33 _ ^22)] 
Ai3 A23 

+ 913923 ^^12^ ^12^(^33 _^t„). (3.28) 
A31 

Note that both a^a and S^^ are now integrals of motion. The two first terms correspond to trivial free atomic 
dynamics. The next two terms represent the standard dynamical Stark shift. Finally, the last term describes an 
effective interaction between levels 1 and 2. The remarkable point is that there is a population transfer (and not only 
phase transfer, as it could be expected from a dispersive interaction) between these two levels without exchange of 
photons. The intensity of the transition 1 <-!■ 2 depends on the difference between the population of level 3 and the 
photon number. Thus, no population transfer between levels 1 and 2 will occur in the sector where the number of 
photons is exactly equal to the initial population of the level 3. It is easy to observe that the (effective) transitions 
1 2 are stronger when A31 = A32; i.e., when the levels 1 and 2 have the same energy (Zeeman-like systems). 
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IV. EFFECTIVE HAMILTONIANS FOR NONLINEAR SU(N) DYNAMICS 
A. Multilevel systems interacting with a quantum field 

Having demonstrated the role played by nonlinear algebras in the systematic construction of effective Hamiltonians, 
we would like to pursue here a natural extension to multilevel systems. 

Specifically, we are interested in considering Hamiltonians that can be represented in terms of the su(N) algebra, 
which naturally arises when describing systems with N relevant levels. Obviously, the evolution of a collection of A 
identical iV-level systems (for definiteness, we assume a cascade configuration, s uch that Ei < Ej for i < j) interacting 
with a single-mode quantum field can be modeled by a Hamiltonian as in Eq. (3.1) with 

■f^ficid = LUfa^a, 

N 

iJatom = (4.1) 

1=1 

H,^t = ^.g,:(a5;*+i+at5!!+i), 

where S" (z = 1, . . . , N) are population operators of the ith energy level, and describe transitions between levels 
i and j. The operators S"-* form the u(N) algebra and satisfy the commutation relations 

By introducing inversion-like operators 



S-J^+i = _ s''), (4.3) 



then (S*^, 5"^*+^) turn out to be the su(N) algebra. 
The Hamiltonian (4.1) admits the integral of motion 



N-l 



N = a^a+J2f^^^Sl^^\ (4.4) 

i=l 

with fii — i{N — i). We also introduce the detunings by 

A, = Ej ~E,- (j - l)c^f , (4.5) 

and we shall assume that Aj satisfy the following resonant condition 

An = 0, (4.6) 

which means that the field in a {N — l)-photon resonance with the atomic system; i.e., E^ ~ Ei — {N — l)uji- Thus, 
the Hamiltonian (4J) can be recast a,s H = Hq + Hint, with 



Ha = UfN + EA 
Hint = ho + V, 



(4.7) 



where E = {En + Ei)/2 and 



N N-1 



A,S", V=J2 9^iaS!l+^ + a^S'l+^). (4.8) 



From our previous experience it seems rather obvious that the operators 

X" = S'\ X'^=aS%, X!f'=at5!_^ (4.9) 
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form a polynomial deformation of su(N). In consequence, and according to our general scheme, we can introduce the 
transformation 



U = exp 



'N-l 



J2 - X"+^] 



where 



Ai+i 



(4.10) 



(4.11) 



will be supposed to be small numbers, <C 1, which means that the transitions are far from the one-photon resonance 
(Ai+i — Ai = Ei^i — Ei — uj{ ^ Qi). Thus, all one-photon transitions are eliminated by ( 1.1C| ) and the transformed 
Hamiltonian takes the form 



N-2 JV-n-1 

En 
(n + l)\ 

Here, the effective interaction constants A" can be obtained from the recurrence relation 



iV— Tl— 1 



A 



(«+i) _ 



'i+nAj 



(4.12) 



(4.13) 



with the initial term A^"* = gi- It is easy to see that A^""''^'' <C A^"''. 

The piece contains only diagonal terms in the atomic operators and depends on the integral of motion N (or, 
equivalently, depends only on the photon- number operator a^a). This operator appears naturally represented as 
an expansion in the small parameter Si whose first term is 



N-l 



(4.14) 



The essential point for our purposes is that, given its structure, this diagonal part cannot be removed from the effective 
Hamiltonian ( [1.12| ). On the contrary, the operator h^d contains only nondiagonal terms that can be eliminated 
by rotations of the type (I-IC) unless some specific resonance conditions are fulfilled. In this respect, let us note 



that the price we pay for eliminating one-photon transitions is the generation of all possible fc-photon transitions 
{k = 2,...,N-l). 



B. Three-photon resonance 



Let us consider the particular case of four-level systems {N = 4) and suppose that there are no transitions in 
one- and two-photon resonance with the field. After eliminating one-photon transitions the transformed Hamiltonian 
(4.12) has the form 



ffW =ho + hd + Kd + ^Af' [a'Sl^ + a^'sl^) + 1 ^ Af ^ [a^S^+' + a^' S^+^ 



(4.15) 



where the interaction constants are defined, according to (4.11) and (4.13), as 

9i _ .92 _ 53 

As' 



ei 



A 



(2) 



51 52 



A2' 
2A2 - A3 

A2(A3-A2)' 



£3 



(4.16) 



A 



(2) 



2A3 



5253 



A3(A2-A3)' 



, (3) ^ 3giff2g3 

' A3A2 ' 



and the resonance condition A4 = (that is, three-photon resonance E4 — Ei — 3uj{) has b een imposed. 

According to the general scheme, the term representing two-photon transitions in (4.15) can be removed using a 
transformation analogous to (4.10): 



U2 — exp 



(2)^„2 c,»+2 _ 



(4.17) 
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where 



J2) _ 



H+2 



(4.18) 



(2) 

is a small parameter because there are no resonant two-photon transitions (Aj_|.2 ~ = £'i+2 — Ei — 2uj{ S> ) and 

thus a^-^^ <C Sj. It is worth noting that the transformation (4.17) does not introduce new terms of order to the 
effective Hamiltonian. 



The diagonal part in (4.15) is 



3 

E 

1=1 



5,e, [ata(S'*+i'+i _ S'') + {S'' + 1)S 



■i+li+l] 



while the nondiagonal term deserves a more careful analysis. Its explicit form is 

hnd 



23 q12\ 



(4.19) 



(4.20) 



It is clear that the first term in the above expression describes a resonant dipole-dipole interaction, under the condition 
A2 = —A3. On the other hand, the second term describes a resonant interaction whenever 2A2 — A3, which is 
incompatible with the pervious one and the absence of one- and two-photon resonances. If no one of these conditions 
are fulfilled, the term h^d can be eliminated by the transformation 



exp 



where 



A. 



H+l 



(4.21) 



(4.22) 



Then, since S"^^ -I- iS^^ -I- S^^ + S"^^ — A and imposing the condition of the absence of initial population in levels 2 and 
3, we obtain the effective Hamiltonian describing three-photon resonant transitions 



H. 



(2) _ gi.92.g3 



eff 



A2A3 ^"'^+ + "^'-^ " + ^/2)Ma(g?/A2 - 53V A3) + 53V A3] + 



(4.23) 



The effective t hree-p hoton Hamiltonian ( 4.23 ) contains a dynamical Stark shift similar to that appearing in the 
two-photon case (3.22). Nevertheless, the essential difference consists in that in the two-photon case the interaction 
term and the Stark shift are of the same order of magnitude, while in the three-photon case the interaction term is 
one order of magnitude less than the Stark shift. This would lead to essential differences in the evolution of some 
observables. 



C. Multimode fields and resonance conditions 



The formalism developed can also be used to treat the interaction of multifrequency fields with atomic systems. 
This type of interaction is much richer because some additional resonance conditions can be satisfied. To this end, 
let us consider a four- level system in a S configuration interacting with two field modes of frequencies tOa andwb, and 



annihilation operators a and 5, respectively. The Hamiltonian describing this interaction is still of the form (3.1) with 

Hficid = LUaa^a + LObb^b, 
4 

i?atom = (4.24) 
i=l 
3 

Hint = l-^9a^a + 9ub)Sl+^ + {ga^a) + 56,&^)5!!+i] , 



4=1 
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where gai and gu ate the corresponding coupling constants. We assume that the following resonance conditions are 
satisfied: 



E4 — El ~ 3ujb, 



El = 2LUb, 



(4.25) 



while all one-photon transitions are out of resonance. In such a case, the Hamiltonian (4.24) may describe, for 
example, a fifth-order process involving the absorption of three photons from one field and the stimulated emission of 
two photons of different frequency. The effective Hamiltonian describing explicitly these tra nsitions can be obtained 
according to the general method. To this end, we note that now the integral of motion (4^) takes now the form 



3 ' 



(4.26) 



i=l 



where /i^ = i(4 — i), and the Hamiltonian can be recast also as in (4.7) with H^^^ = ^0 + + 14: where 

4 



(4.27) 



3 

E 



ii+l 



X 



Here the polynomial deformation at hand is defined by 



Y^3 



(4.28) 



and S = LUf, ~ uja will be taken, for definiteness, as positive. 

Now we can eliminate all the one-photon nonreso nant t ransitions generated by means of a couple of transformations 
(one for mode a the other for mode b) identical to (4.10). The transformed Hamiltonian (up to order 1/A^) is 



H. 



(1) 



(4.29) 



where h^^^ and h''^' are defined according to (1.14) for both modes and the coupling constants XY' and A^'^'' are given 
by ( pT6| ). The constant is 

.9a3562 gb3ga2 



,(3) 



Aab) 

A4 - A3 A3 - A2 ■ 

This effective Hamiltonian deserves some comments. First of all, when the resonance condition 

E4 — E2 = + CiJf, 



(4.30) 



(4.31) 



is fulfilled (which is compatible with the two- and three- photon resonance conditions), then the last term in the 
Hamiltonian describes resonant transitions between levels 2 and 4, as a result of the simultaneous absorption and 
emission of quanta from both modes, and take place only if levels 2 and 4 are populated. In absence of such an 

additional resonance condition, the Hamiltonian (4.29) describes two simultaneous and competing processes: a first- 

(2) 

order process {X\ ~ V^) involving two-photon transitions between levels 1 and 3, and a second-order process 
(A^'^'' ~ 1/A^) of three-photon transitions between levels 1 and 4. 



FURTHER EXTENSIONS AND CONCLUDING REMARKS 



We have developed our method putting special emphasis in polynomial deformations of algebras su(7V), due to 
the outstanding role that they play in describing the interaction of A^-level systems with quantum fields. However, 
by no means the method is re strict ed to such kind of algebras. In fact, many other phenomena can be modeled by 
Hamiltonians quite similar to (2.11), namely 



a Yo + g{Y++Y^) + C, 



(5.1) 
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where C is some integral of motion and a a constant. In these theories, the polynomial deformation is defined in the 
following fashion in the Cartan-Weyl basis: 



where $(10) a-re appropriate structure polynomials. A detailed study of the applications of these algebras for solving 
evolution problems in nonlinear quantum models may be found in Ref . | [T9[ | . For these model, the machinery of small 
rotations works well and constitutes the most systematic ways of constructing effective Hamiltonians. 

In summary, what we expect to have accomplished in this paper is to develop an appropriate and systematic tool 
for obtaining effective Hamiltonians that describe nonlinear optical phenomena. Our approach exploits the existence 
of deformed algebras (arising as dynamical symmetries of the corresponding process) to construct small nonlinear 
rotations that perform the task. 

This method constitutes a firm algebraic setting free from some inconsistencies found in the traditional adiabatic 
elimination of variables and provides unambiguous recipes for proceeding with any model Hamiltonian. 
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